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CURVILINEAR SURFACE* 


By V. F. Droblenkov 


A number of semi empirical approximate methods exist for determining 
the characteristics of the turbulent boundary layer on a curvilinear 
surface. At present, among these methods, the one proposed by L. G. 
Loitsianskii (ref. 1) is given frequent practical application. This 
method is sufficiently effective and permits. In the case of wing pro- 
files with technically smooth surfaces, calculating the basic character- 
istics of the boundary layer and the values of the overall drags with 
an accuracy which suffices for practical purposes. 

The idea of making use of the basic integral momentum equation 


+ Y (2 + h) = ISL 
a* V pV 2 


( 1 ) 


proves to be fruitful also for the solution of the problems in the 
determination of the characteristics of the turbulent boundary layer 
on a rough surface. In fact, using, for instance, the well-known values 
of the averaged longitudinal velocities of the turbulent layer on a 
plate with a developed roughness 


22. = 2.5 In -2- + 8.48 

v * k s 


the equation for the relationship between the local-friction coeffi- 
cient Cf and & KX 'I k B may be obtained in the form 



*"Turbulentnyi pogranichnyl slot na sherpokhovatoi krivolineinoi 
poverkhnosti." Izvestiia Akademii Nauk SSSR, Otdelenie tekhnicheskikh 
nauk [Reports of the Academy of Sciences of the S.S.S.R., Department of 
Technical Sciences^], no. 8, 1955; PP* 17-21* 



2 


NA.CA. m 1440 


Here and above , 


the following 


symbols are used: 


c f 



8^ = 


6 X = 



H = 


Sf 

V 


v* 



where S 305 " is the momentum thickness loss, 6 X the displacement thick- 
ness, k s the height of the protuberances of the equivalent roughness, 

Vq the velocity of the approach flow, V the velocity on the upper 

limit of the boundary layer, and v® the velocity of the inside layer. 

Figure 1 represents L]_ = lg ^pVq 2 / Tq^ as a function of = lg R 301 

and = lg^S 30 ^ k s ^) where curve 1 pertains to the laminar layer and 
curve 2 to the turbulent boundary layer for which 


-I2_ , o.oo655(b xx )' 1/ ' 6 

pV 


Both curves refer to the case of a technically smooth surface, whereas 
curve 3 applies to the turbulent boundary layer 
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T o 

pv 0 £ 


0.003 



-1/6 


(3) 


in the case of a rough surface. 

As can "be seen from figure 1, the equality (3); replacing the 
equality (2), is correct over the wide range of variation of the rela- 
tive roughness^- and is more convenient for practical calculations. With 
the aid of the equality (3)* the solution of the fundamental houndary- 
layer equation (l) no longer presents any difficulty for rough plane 
surfaces or for rough curvilinear surfaces. 

Using the equality (3), the fundamental momentum equation (l) for 
the case of a longitudinal flow about the plates (V = V v ’ - °) 

as** _ 

'pv 0 2 


is rewritten in the form 


dS** 

dx 


0.0031 



-l/6 


By integration the following equation is obtained: 



- 1 0.003l(k s ) 


1/6 

x + C 


Taking into consideration that the turbulent boundary layer is 
established from the leading edge ^Xq = 0, Sq^ = 0^, C = 0 is 

obtained, and finally 


-krf necessary, the equality (2) can be approximated by an expression 
which is, as to form, analogous to the equality (3), but reflects much 
more strongly the dependence of (2) in the case of essentially different 
values of k s and Cf. 
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6* ** = 0.0080(k s ) l/7 - 6/7 


x ' ’ or = 0.008o(^~ 


,- 1/7 


V k s 


According to the value 5 , the local turbulent-friction coeffi- 

cient on an entirely rough plate can he found as 


c f 




and the value of the overall drag coefficient of a rough plate of a 
given length L as - 


Cf = 


■Jo T ° - So 0f <5) ■ °- 0l62 fe)" 1/7 (4) 


i p v o 2l 


i pV 


The final results of the calculations, carried out according to 
the equality (4), are plotted in figure 2 where the solid line repre- 
sents Cf as a function of lg(L/k s ), obtained by calculation according 
to formula (4), and the dotted line represents Schlichting's data. As 
can be seen, the calculated results show good agreement with the proved 
values of the overall turbulent-friction coefficients of plates for the 
regime of developed roughness obtained by Schlichting. This latter 
fact permits us to consider the obtained solutions as correct over a 
wide range of variation of relative roughness. 

With the aid of equality (5) and the momentum equation (l) the 
turbulent boundary layer can also be calculated on a rough surface in 
the case of a longitudinal pressure drop. T 

Let it be assumed that for the case of a turbulent boundary layer 
on a rough surface, there exists the function Gsh(^sj S 30 ^ which has 


*TTACA reviewer^ s note: The subscript sh stands for the Russian 

word "sherpokhovatoi, " which translated, means "rough." 
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the property that, when multiplied hy the quantities V' and 

T 0 /pV 2 , it "becomes independent of the Reynolds number and the parameter 


f = 


V'&xx 


r sh 


( 5 ) 


and its function 

W f ) = \ G sh = £sh (6) 

pV^ 

/yy 

5 

Determination of the form of the function Gg^ is carried out here 

also, as in L. G. Loitsianskii' s method, "by starting from the analogy 
with the 1 nmlnsr boundary layer, under the assumption that for all 
values f the form of the function G sll coincides with that for a 

turbulent flow, in the case of a longitudinal flow about a rough plate. 

Since V' is equal to 0 on the plate, according to the equali- 
ties (5) and (6) on the plate, f = 0 and = const. Then, assuming 

for simplification (; sll = const. = 1, according to the equalities (6) 

and (5) 


1 = 


‘0 


pv 0 


7J- Ggk = 0.0031^g— - 


- 1/6 


G sh 


Hence, there is found, also, the form of the function G g -^ : 

- ■ 

The value of the argument of the function G B h for every section 
of the boun dar y layer of the surface under consideration is determined 
from the solution of the fundamental equation, (l) . This solution may 
be represented by u an expression analogous to the solution obtained by 
L. G. Loitsianskii for a turbulent regime of flow about a technically 
smooth surface . Following his method, after several transformations, 
instead of the equality (l), an ordinary differential equation can be 
obtained for determining f on the rough surface 
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df V 


dx 


= JSLL F(f ) + f ^ 

V 


(8) 


where 


F(f) = (l + m)£ sh - f[5 + an+(l + m 




Since the solution of equation (8) is known, it"is not difficult 
to find 5 X , said t o/p^O^ ^ ^ orm °f the function Gsh for 

every section of the turbulent boundary layer on a rough surface is 
known, as well as that for the function F(f). 

If the longitudinal drop is small, the values of f also will be 
small, and it will not involve a great inaccuracy to substitute for H 
and £ s k, in the equality for F(f ), their values for f' = 0. 

Then, assuming that H(f) Is equal to 1.3 for f = 0 (for large 

Reynolds numbers and x/kg > 1 x 10 / and £ s k Is equal to 1.0, it 
follows that 

F(f) « a - bf = 1.17 - k.85£ 


and the solution of the fundamental equation of the turbulent boundary 
layer on a rough surface is used for determining the quantity f(x) 
in the form of a simple quadrature ~ 
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or 


f<x) ' a &C * wfefo v5 ' 8?to 

where the turbulent boundary layer begins directly at the leading edge. 

The calculation of the values of f(x) according to the values of 
the velocities V on the upper limit of the boundary layer does not 
present any difficulties. 

With the aid of the equalities ( 5 ) and (j) the values of the 
momentum thickness loss at the trailing edge are found to be 


8k** 



V b_1 dx 


1 

m+1 


Consequently, the value of the total-drag coefficient of a body with a 
rough surface in a two-dimensional, completely steady, turbulent flow 
regime becomes 



For the purpose of additional verification and proof of the correct- 
ness of the resulting conclusions regarding the equalities ( 3 ), (4), 
and (9), the values C x of the NACA 0012 wing profile were calculated 
for three surface conditions: technically smooth, rough with the value 

L/ks = 8.12 x 10^, and rough with Lfk s = 2.03 X 10^. The results of 
the calculation were compared with the results of tests in a variable- 
density tunnel with a wing of that profile, an aspect ratio of 6 , and a 
chord length of 203 mm. The test data and calculated results are pre- 
sented in figure 3* There, the curve X corresponds to the drag coeffi- 
cient of a technically smooth plate in the case of a completely turbulent 
boundary layers to the measured values of C x for the wing profile, 
curve II corresponds to an accurately handmade surface structure; and 
curves III and IV represent conditions for surfaces covered with grains 
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of carborundum powder of the size k g = 0.025 nm and' k g = 0.01 mm, 

respectively. As can be seen from the data of figure 5, the calculated 
results (solid lines 1, 2, and 3) and the test results (dotted lines) 
agree with each other sufficiently well for high Reynolds numbers when 
a turbulent regime of flow is established along practically the entire 
profile length. The insignificantly excessive values of the calculated 
quantities C x (obtained for a completely turbulent flow regime) com- 
pared to the test quantities is most pronounced in the case of a tech- 
nically smooth surface condition of the profile, and it may be assumed 
that it is the result of the influence of the laminar flow regime which 
occurs in test6 in the region of the leading edge. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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